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A  new  relationship  between  the  roughness  height  and  the  main  hydrodynamic  and  sediment  parameters 
for  plane  beds  under  oscillatory  conditions  is  presented.  In  order  to  derive  such  a  relationship,  a  large 
data  base  encompassing  plane-bed  experiments  was  compiled  from  previous  investigations  and 
analyzed.  Different  methods  to  estimate  the  roughness  height  were  investigated.  Comparisons  between 
the  data  and  different  existing  predictive  formulas  for  the  bed  roughness  obtained  from  the  literature 
were  also  performed.  A  relationship  involving  the  grain  size  and  Shields  parameter  only,  which  is 
commonly  proposed,  appeared  to  be  insufficient  to  characterize  the  roughness.  The  roughness  height 
was  also  found  to  be  a  function  of  the  sediment  density  and  the  settling  velocity.  A  critical  Shields 
parameter  was  identified  up  to  which  the  effective  roughness  ratio  is  proportional  to  the  Shields 
parameter  only.  The  new  empirical  equation  developed  in  this  study  provides  the  highest  predictive  skill 
for  all  conditions  investigated. 

©  2008  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Accurate  sediment  transport  rate  predictions  in  coastal  areas 
closely  depend  on  estimates  of  the  resistance  to  oscillating  flows  on 
the  sea  bottom.  A  roughness  height  is  typically  employed  to  char¬ 
acterize  the  resistance  properties  of  the  bottom.  Although  the 
roughness  height  remains  quite  difficult  to  determine,  it  is 
a  fundamental  parameter  for  calculating  sediment  transport  rates. 
Numerical  models  to  compute  sediment  transport  rate  in  river  or 
coastal  environments  often  employ  relationships  which  depend  on 
the  Shields  parameter.  The  bed  roughness  is  an  important  input 
quantity  in  the  Shields  parameter,  and  it  is  usually  unknown  even 
for  the  simpler  case  of  steady  current  conditions. 

In  general,  for  dimensional  reasons,  the  bed  roughness  height  of 
a  flat  and  fixed  bed  is  given  in  terms  of  the  Nikuradse  roughness 
height  (ks).  For  flat  beds,  it  is  expected  to  be  on  the  order  of  the 
median  grain  diameter  or  of  some  larger  grain  size  percentiles 
(/<s  =  l-5  times  dso,  (fe,  ds4,  or  dgo  according  to  the  literature). 
However,  the  value  of  ks  varies  considerably  depending  on  the 
configuration  of  the  grains  forming  the  roughness  of  the  flow 
boundaries.  At  high  Shields  parameter  values,  the  sediment  moves 
along  the  bottom  mainly  in  a  layer  denoted  as  sheet-flow.  The 
presence  of  a  relatively  thin  sheet-flow  layer  with  high  sediment 
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concentration  markedly  affects  the  flow  above  it.  One  important 
aspect  of  the  presence  of  a  sheet-flow  layer  is  the  increased 
roughness,  compared  to  a  situation  without  such  a  layer.  Under 
sheet-flow  conditions,  the  roughness  height  may  be  several  orders 
of  magnitude  larger  than  for  a  fixed  bed.  This  is  probably  caused  by 
the  increased  energy  dissipation  in  the  sheet-flow  layer  due  to 
interaction  between  individual  sediment  grains  as  well  as  between 
the  sediment  and  the  fluid.  It  is  often  assumed  that  the  roughness 
height  under  sheet-flow  conditions  is  on  the  order  of  the  layer 
thickness  (Wilson,  1987;  Van  Rijn,  1993). 

Wilson  (1966)  performed  a  series  of  tests  in  a  pressurized 
conduit  in  order  to  obtain  high  shear  stress  values.  Based  on  these 
results,  he  proposed  a  linear  relationship  between  the  relative 
roughness  and  the  Shields  parameter.  Some  more  recent  studies 
(Nnadi  and  Wilson,  1992;  Sumer  et  al.,  1996)  confirm  this  rela¬ 
tionship  with  the  Shields  parameter.  However,  all  these  experi¬ 
ments  were  performed  for  steady  currents.  Few  experiments  have 
been  performed  under  oscillatory  flow  conditions,  and  mostly  in 
the  laboratory.  Moreover,  if  the  shear  stress  at  the  bottom  (and  thus 
the  roughness  height)  can  be  easily  estimated  in  the  case  of  steady 
flows  using  energy  slope  measurements  (Camenen  et  al.,  2006), 
this  estimation  becomes  more  difficult  in  the  case  of  oscillatory 
flows.  An  extensive  laboratory  data  set  was  compiled  in  this  study, 
allowing  for  a  more  comprehensive  analysis  of  this  difficulty.  Five 
methods  were  compared  to  estimate  the  roughness  height  for 
oscillatory  flows  based  on  the  velocity  profiles,  energy  dissipation, 
momentum  equation  in  the  fluid/sediment  mixture,  sheet-flow 
layer  thickness,  and  bed-load  sediment  transport.  Most  of  the 
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Nomenclature 

z0 

roughness  lenght  (z0  =  fcs/30), 

Oib ,  0is 

coefficients, 

wave  orbital  semi-excursion  at  the  bottom, 

Pb 

coefficient, 

C 

sediment  concentration, 

5 

thickness  of  the  viscous  (Stokes)  layer  5  =  y/ vTw/n , 

d 

sand  diameter  (grain  size), 

db 

thickness  of  the  erosion  depth, 

^50 

median  sand  diameter, 

ds 

thickness  of  the  sheet-flow  layer, 

d* 

dimensionless  grain  size, 

K 

von  Karman  constant, 

De 

time-averaged  rate  of  energy  dissipation  due  to  bed 

V 

kinematic  viscosity  of  the  fluid  (v  =  10  6  m2/s  for  clear 

friction, 

water), 

fe 

energy  dissipation  factor, 

Vj 

eddy  viscosity, 

fw 

wave  friction  coefficient, 

(D 

wave  angular  frequency  (w  =  2tc/Tw), 

g 

acceleration  due  to  gravity, 

(pr 

critical  dynamic  angle  of  internal  friction  of  sediments, 

ks 

roughness  height, 

P 

water  density, 

p 

pressure, 

Pm 

fluid-sediment  mixture  density, 

Ppl 

phase-lag  parameter, 

Ps 

sediment  density, 

rw 

wave  asymmetry, 

Tw 

bed  shear  stress  due  to  wave  oscillatory  flow, 

t 

time, 

6 

Shields  parameter  due  to  wave  oscillatory  flow, 

Tw 

wave  period, 

Merest 

maximum  Shields  parameter  due  to  wave  oscillatory 

u 

horizontal  velocity, 

flow  at  the  crest  of  the  wave, 

uw 

wave  horizontal  velocity  at  the  bottom, 

^trough 

maximum  Shields  parameter  due  to  wave  oscillatory 

u* 

friction  velocity, 

flow  at  the  trough  of  the  wave, 

li  00 

free-stream  velocity, 

S' 

skin  Shields  parameter  due  to  wave  oscillatory  flow, 

U 

velocity  defect  function, 

6 

mean  value  of  the  the  absolute  value  of  the  Shields 

Uw 

wave  orbital  velocity  at  the  bottom, 

parameter  due  to  wave  oscillatory  flow  over  a  wave 

Uw,  cr 

critical  wave  orbital  velocity  at  the  bottom  for  the 

half-period, 

inception  of  sheet-flow, 

Ocr 

critical  Shields  parameter  for  the  inception  of  motion, 

X 

horizontal  coordinate, 

@cr,  ur 

critical  Shields  parameter  for  the  inception  of  upper- 

W 

vertical  velocity, 

regime, 

Ws 

settling  velocity, 

@cr,  sf 

critical  Shields  parameter  for  the  inception  of  sheet- 

z 

vertical  coordinate, 

flow, 

existing  methods  are  actually  indirect  as  the  bed  shear  stress  is  not 
directly  measured  but  estimated  from  other  parameters.  Thus, 
some  additional  uncertainties  exist  based  on  the  theoretical 
formulation  to  calculate  the  bed  shear  stress  from  these  other 
parameters.  Bayram  et  al.  (2001)  encountered  the  problem  of 
finding  reliable  relationships  to  predict  roughness  height  during 
sheet-flow  conditions  in  the  surf  zone  when  evaluating  the 
predictive  capability  of  available  longshore  sediment  transport 
formulas  against  field  data,  which  motivated  the  present  study. 

The  term  “equivalent”  roughness  height  (or  bed  shear  stress)  is 
used  here  as  it  does  not  necessarily  correspond  to  the  real  value  but  to 
the  one  obtained  in  a  “clear  water”  model  over  a  fixed  bed  to  simulate 
the  same  flow  discharge  and  hydraulic  gradient  as  measured  for  the 
mobile  bed  situation.  The  analogy  with  an  oscillating  turbulent  clear 
water  flow  over  a  fixed  bed  is  assumed  in  order  to  use  classical 
formulas  for  the  friction  coefficient  such  as  the  Jonsson  (1966) 
formula.  By  definition,  the  equivalent  bed  shear  stress  is  thus 
different  as  the  real  bed  shear  stress  below  a  sheet-flow  layer. 

Nielsen  (1992)  investigated  the  equivalent  Nikuradse  roughness 
of  natural  sand  beds  exposed  to  both  steady  and  oscillatory  flows. 
Based  on  data  collected  by  Carstens  et  al.  (1969)  and  Lofquist  (1986), 
Nielsen  showed  that  flat  beds  in  oscillatory  sheet-flow  (z.e.,  Shields 
parameter  0  >  0.8)  generally  exhibit  roughness  values  on  the  order 
of  100-200  grain  diameters.  However,  steady  sheet-flow  data  by 
Guy  et  al.  (1966)  produced  roughness  heights  one  order  of  magni¬ 
tude  smaller  than  the  oscillatory  flow.  No  satisfactory  explanation 
has  been  given  for  this  difference  between  the  roughness  values  of 
sand  beds  under  steady  and  oscillatory  flows.  Nielsen  (1992) 
reasoned  that  the  observed  additional  energy  dissipation  under 
waves  due  to  percolation  might  explain  the  difference. 

The  objectives  of  the  present  study  are  twofold:  (1)  to  comple¬ 
ment  the  findings  of  Nielsen  (1992)  by  analyzing  additional 


high-quality  available  data  sets;  and  (2)  to  propose  a  new  predictive 
equivalent  roughness  height  relationship  for  practical  applications 
in  sediment  transport  modeling  in  coastal  areas,  using  various  data 
sets  from  the  literature.  In  developing  an  improved  predictive  rela¬ 
tionship  for  the  roughness  height,  several  existing  formulas  were 
also  compared  with  the  compiled  data  set.  The  first  part  of  the  paper 
presents  the  previous  studies  (c/.  Section  2),  the  different  methods  to 
estimate  the  roughness  height  (c/.  Section  3.1 ),  the  data  sets  used  in 
the  present  study  (cf.  Section  3.2),  and  the  validation  of  the  methods 
to  estimate  the  roughness  height  (cf.  Section  3.3).  The  effects  of  the 
main  parameters  on  the  roughness  height  and  the  new  predictive 
relationship  are  discussed  in  the  second  part  (cf  Section  4). 

2.  Roughness  height  relationships 

In  general,  the  roughness  height  is  assumed  not  to  vary  with 
time  even  for  oscillating  flows.  Thus,  all  existing  formulas  yield  an 
average  value  of  the  equivalent  roughness  height.  Recently,  some 
authors  (Zala  Flores  and  Sleath,  1998;  Nielsen,  2002)  showed  that  it 
may  vary  around  its  mean  value.  However,  the  purpose  of  the 
present  paper  remains  to  improve  the  prediction  of  its  mean  value. 

2.1.  Existing  relationships 

In  general,  the  main  problem  in  bed  roughness  prediction  under 
moving  bed  conditions  is  that  the  roughness  height  depends  on  the 
flow  variables  as  well  as  on  the  sediment  transport  rate.  In  the 
literature,  several  existing  relationships  for  equivalent  roughness 
under  sheet-flow  conditions  have  been  proposed  (see  Table  1). 
Most  of  the  existing  formulas  assume  a  relationship  between  an 
effective  roughness  ratio  (ratio  between  the  effective  roughness 
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Table  1 

Relationships  to  predict  equivalent  roughness  for  oscillatory  sheet  flow  as  they 
chronologically  appeared  in  the  literature. 


Grant  and  Madsen  (1982) 

Wilson  (1989) 

Nielsen  (1992,  pp.154-155) 
Madsen  (1993) 

Xu  and  Wright  (1993) 
Ribberink  (1998) 


k s_ 

^50 


_k< 

^50 


k s_ 

^50 

k s_ 

^50 

k s_ 

^50 

k s_ 

^50 


=  50 


160(s  +  0.5)6cr  j-  —  0.7J  (2) 

(3) 

7oVd 


=  15 

-  15  d 

max[3, 1  +  6(<  8  >  -1)] 


(4) 

(5) 

(6) 
(7) 


height  ks  and  a  characteristic  grain  size  d)  and  the  Shields  param¬ 
eter  based  on  the  maximum  shear  stress  6 ,  according  to, 


on  the  instantaneous  shear  stress.  Their  measurements  show  some 
influence  of  pressure  gradient  and  inertia  forces  for  large  shear 
stresses  where  a  type  of  plug  flow  is  observed.  They  defined  the 
parameter  S  representing  the  ratio  of  inertia  to  gravity  forces  acting 
on  individual  grains  to  characterize  this  phenomena: 


_  Uw(o 

~  (sfl)g 


(11) 


From  the  authors’  point  of  view,  the  Shields  parameter  better 
characterizes  the  ratio  of  inertia  to  gravity  forces.  Also,  the  model 
proposed  by  Zala  Flores  &  Sleath  links  the  ratio  <y(Aw/w)  to  S  which 
directly  leads  to  the  Shields  parameter  for  a  linear  case.  For  S  >  0.2, 
they  observed  a  phase  lead  of  the  fundamental  component  of  the 
velocity  which  varies  depending  on  the  relative  rates  of  erosion  and 
deposition  of  the  sediments.  This  confirms  the  importance  of  the 
settling  velocity  and  density  of  the  sediments  on  the  roughness 
height.  Their  measurements  showed  that  the  velocity  profiles  are 
different  during  the  acceleration  phase  and  the  deceleration  phase, 
which  may  induce  a  net  sediment  transport. 

Nielsen  (2002)  showed  that  the  acceleration  asymmetry 
enhances  the  bed  shear  stress.  Fie  proposed  to  use  a  model  for  the 
skin  shear  stress  that  generates  stronger  shear  stresses  under  the 
more  abruptly  accelerated  part  of  the  instantaneous  free-stream 
velocity, 


where  0  varies  between  0.5  and  2  depending  on  the  authors. 

In  Table  1  d^o  is  the  median  diameter,  and  O'  and  6  are  the  skin 
and  total  Shields  parameters,  respectively.  The  Shields  parameter  is 
defined  by, 


in  which  u*  is  the  bed  friction  velocity,  s  is  the  relative  density  of 
sediment  grains  (s  =  ps/p,  where  ps  and  p  are  the  density  of  sedi¬ 
ment  grains  and  fluid  (water),  respectively)  and  g  is  the  accelera¬ 
tion  due  to  gravity. 

The  shear  velocity  can  be  calculated  with  the  aid  of  the  friction 
factor  defined  by  Jonsson  (1966), 


Mt)  =  (s-iterf50[f  (M0,  vt)\f  (M0,  <pt) 

f(Uw(t),(pt)  =  cos<ptuw(t)  +  smcpt — 1 ^ - - 

where  the  friction  coefficient  fw  is  calculated  using  the  skin 
roughness  height  (/<s  ~  2dso),  (pt  is  the  phase  shift  between  free- 
stream  velocity  and  bed  shear  stress  at  the  peak  frequency  (Nielsen 
found  cpt  ~  37.5°),  and  5t  is  the  time  step  (c/.  Eq.  (7)  with  n  =  0  in 
Nielsen,  2002,  where  ne  [0;  1  ]  is  an  interpolation  parameter). 

Camenen  and  Larson  (2007a)  proposed  a  simple  model  which 
increases  the  friction  coefficient  (previously  assumed  constant  over 
a  wave  period)  for  accelerated  flow  and  decreases  it  for  decelerated 
flow.  The  mean  Shields  parameters  in  the  onshore  (0w(t)  >  0)  or 
offshore  direction  ( 6w{t )  <  0)  are  thus  modified  as  follows, 


214  Tw  r 

/W  U2  1/2  p(AwM)2  W 

where  Uw  is  the  maximum  bottom  wave  orbital  velocity,  tw  the 
wave-induced  bed  shear  stress,  Aw  the  amplitude  of  the  near-bed 
wave  orbital  motion,  and  (o  =  2n/Tw  the  angular  frequency  (Tw  is  the 
wave  period).  An  assumption  usually  made  for  the  coastal  envi¬ 
ronment  is  that  the  rough  turbulent  regime  is  fully  developed. 
Thus,  the  friction  factor  could  be  estimated  by  the  formula  sug¬ 
gested  by  Swart  (1974),  which  is  an  explicit  approximation  to  the 
implicit,  semi-empirical  formula  given  by  Jonsson  (1966): 


fw  =  exp 


5.2  -r^ 


fc'0-194 


-5.98 


(10) 


^w, on*/off*  —  ^w,on/off0  ^  ra)  (13) 

with  ra  =  (Ja  -  Td)/Td  being  the  acceleration  asymmetry  (Ta  and  Td 
are  the  half  periods  where  the  acceleration  is  positive  and  negative, 
respectively).  The  Camenen  and  Larson  model  leads  to  very  similar 
results  compared  to  the  Nielsen  equation  except  when  ra  >  0.7 
where  the  Nielsen  model  tends  to  diverge.  Therefore,  the  acceler¬ 
ation  may  be  assumed  to  significantly  affect  the  instantaneous 
friction  coefficient  and  roughness  height  (and  thus  the  net  sedi¬ 
ment  transport)  but  not  their  mean  values.  The  equation  proposed 
in  this  study  (i.e.  Eq.  (33))  should  thus  be  valid  independently  of  the 
acceleration. 

3.  Oscillatory  flow  data  set 


Jonsson  (1980)  suggested  an  upper  limit  of  0.3  for  the  value  of/w. 
The  skin  shear  stress  is  computed  following  the  presented  method 
and  assuming  that  ks  =  2dso.  The  choice  of  ks  =  2dso  is  based  on 
a  previous  study  for  steady  flows  (Camenen  et  al.,  2006).  It  corre¬ 
sponds  to  a  mean  value  for  the  lower-regime  where  the  sediment 
transport  does  not  affect  the  roughness  height. 


In  order  to  investigate  roughness  height  under  plane-bed 
conditions,  a  wide  range  of  existing  data  sets  were  compiled  and 
analyzed.  As  previously  mentioned,  a  major  difficulty  in  the  case  of 
oscillatory  motion  is  to  determine  the  bottom  shear  stress.  Several 
methods  to  estimate  the  total  bottom  shear  stress  are  presented 
following  a  literature  review. 


2.2.  Acceleration  effects  on  the  roughness  height 

Zala  Flores  and  Sleath  (1998)  were  the  first  who  observed 
acceleration  effects  on  the  instantaneous  velocity  profile,  and  thus 


3.1.  Methods  to  estimate  bed  roughness  and  total  shear  stress 

Since  at  the  moment  there  are  no  possibility  to  directly  measure 
the  shear  stress  over  a  moving  bed,  several  methods  are 
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investigated  to  estimate  the  bed  shear  stress.  As  observed  before, 
these  methods  lead  to  the  estimation  of  the  equivalent  bed  shear 
stress  in  “clear  water”  to  correctly  reproduce  the  flow  over  a  bed 
assumed  fixed  and  to  compensate  for  the  differences  with  the  real 
two-phase  flow. 


3.1.1.  Velocity  profile 

The  velocity  profile  in  the  wave  boundary  layer  can  be  computed 
using  the  turbulent  boundary  layer  momentum  equation  for 
unsteady  flow.  Assuming  that  the  flow  inside  the  boundary  layer  is 
essentially  horizontal  (w  =  0),  the  unsteady  boundary  layer 
approximation  to  the  Reynolds-averaged  horizontal  momentum 
equation  with  an  eddy  viscosity  closure  has  the  form, 


du  du 

di  +  U  to 


1  dp 
p  dx 


(14) 


where  x  and  z  are  the  horizontal  and  vertical  coordinates,  respec¬ 
tively,  u  is  the  horizontal  velocity,  t  the  time,  p  the  pressure,  and 
vj (z)  the  eddy  viscosity.  This  can  be  coupled  to  a  flow  approxima¬ 
tion  for  the  velocity  above  the  boundary  layer  or  free-stream 
velocity  u 


du oo  _  _^dp  (  . 

~dt~  ~  ~~p  6x  1  ] 

Neglecting  the  non-linear  convective  velocity  term  ( udu/dx  ~  0) 
and  recognizing  that  9Uoo/0x  =  O,  Eqs.  (14)  and  (15)  reduce  to  an 
equation  in  terms  of  velocity  difference: 


^(u“-u)  =  i(^(z4(u»-u))  (16) 

If  the  wave  motion  is  assumed  to  be  sinusoidal  (i.e., 
Uoo  -  u  =  Re(Uewt)),  the  flow  velocity  within  the  boundary  layer  can 
be  expressed  by  an  ordinary  differential  equation: 


dz2  vT  dz  vT 

An  analytical  solution  to  Eq.  (17)  (with  the  boundary  conditions 
=  atz  =  Zo  =  ksl 30,  and  U  =  0  for  z» <5W,  where  dw  is  the  wave 
boundary  layer  height)  can  be  obtained  close  to  the  bottom  for 
small  values  of  Zoco/(ku *)  where  the  logarithmic  velocity  law  can  be 
shown  to  be  valid  and  the  eddy  viscosity  is  a  linear  function  of  the 
vertical  position,  vt=ku*z  (Grant  and  Madsen,  1986), 


assuming  that  other  dissipating  mechanisms  than  friction  were 
negligible  in  the  experiments,  the  energy  dissipation  factor  fe  can 
be  predicted  using  the  Jonsson  (1966)  relationship, 


fe  =  _ * _ ,  =  Tw(t)U“(t)  ,  (20) 

2/(3ir)j o(Awoj)3  2/(3n)p(Aww)3 

where  De  is  the  time-averaged  rate  of  energy  dissipation  due  to  bed 
friction,  in  which  rw(t)  is  the  instantaneous  wave-induced  bed 
shear  stress. 

Although  the  two  friction-related  coefficients  fw  and  fe  are 
different  according  to  their  definitions,  Nielsen  (1992,  pp.27-28) 
showed  that  the  experimental  scatter  of  the  measurements  of  one 
or  the  other  over  natural  sand  beds  is  so  large  that  for  practical 
purposes  fw  and  fe  can  be  assumed  equal.  The  roughness  height  can 
be  calculated  inversely  using  Eq.  (10)  once  fe  is  determined  from 
experimental  data. 

3.1.3.  Momentum  equation  in  the  fluid/ sediment  mixture 

The  bed  shear  stress  can  also  be  obtained  from  the  momentum 
equation  for  the  fluid/sediment  mixture.  In  case  of  pure  oscillatory 
flow  with  no  superimposed  current,  the  bed  shear  stress  is  given  by 
Dick  and  Sleath  (1991), 

Tw  =  J  ^(PmU  ~  PU°°)dZ  (21) 

where  pm  is  the  density  of  the  fluid/sediment  mixture, 


pm  =  (l-C)p  +  Cft  (22) 

in  which  C  is  the  volumetric  sediment  concentration. 

However,  prediction  of  bottom  shear  stress  using  the  above 
relation  requires  additional  assumptions  on  the  concentration  C, 
which  varies  in  the  moving  layer  from  0.1  to  0.6  (Wijetunge  and 
Sleath,  1998).  Another  difficulty  is  to  define  the  height  where  the 
bed  shear  stress  tw  is  calculated.  Dick  and  Sleath  (1991 )  proposed  to 
use  the  initial  bed  level.  It  seems  that  the  still  bed  level  during  grain 
motion  is  more  appropriate  but  it  is  not  easy  to  define  (where  the 
maximum  concentration  is  reached).  The  roughness  height  is 
calculated  as  in  Section  3.1.2  where  fw  is  deduced  from  the  defini¬ 
tion  of  the  shear  stress  (c/.  Eq.  (9)). 

This  method  differs  from  the  others  as  it  is  supposed  to  give  the 
“real”  bed  shear  stress  under  the  sheet-flow  layer  and  not  the 
“equivalent”  bed  shear  stress. 


U(z,t)  2  .  ,  x.  z 

;  =  —  simpcos  (cot  +  (p)\n—  (18) 

Uw  TC  Zq 

where  Uw  is  the  maximum  wave  orbital  velocity  and  <p  is  the  phase 
lead  between  free  surface  and  near-bottom  wave  orbital  velocity 
given  by: 


tan  (p 


iz/2 

In  1.15 


(19) 


Thus,  it  is  possible  to  fit  Eq.  (18)  to  an  experimental  data  using  ks 
and  u *  as  they  are  the  only  unknowns.  Since  this  method  corre¬ 
sponds  to  the  definition  of  the  roughness  height  introduced  for 
steady  currents,  it  will  be  subsequently  denoted  as  the  “reference 
method”. 


3.1.2.  Energy  dissipation  factor 

Nielsen  (1992,  pp.146-153)  used  the  experimental  data  by 
Carstens  et  al.  (1969)  and  Lofquist  (1986)  to  estimate  the  wave 
energy  dissipation  factor.  In  these  two  data  sets  are  the  shear  stress 
time  series  (Lofquist)  and  dissipation  measurements  (Carstens 
et  al.)  available  for  sand  beds  exposed  to  oscillatory  flows.  By 


3.1.4.  Sheet-flow  layer  thickness 

In  a  series  of  experimental  studies  with  steady  flow,  Wilson  and 
Pugh  (1989)  and  Pugh  and  Wilson  (1999)  showed  that  the  rough¬ 
ness  height  is  proportional  to  the  thickness  of  the  moving  layer  ds 
(defined  as  the  difference  between  the  top  of  the  moving  bed  level 
and  the  still  bed  level)  with  a  factor  0.5.  Wilson  (1989)  made  an 
analysis  of  the  sheet-flow  friction  during  oscillatory  motion  and 
proposed  a  similar  relationship,  which  may  be  written, 

ks  =  as5s  =  a8d50  9  (23) 

where  as  =  0.5  and  ao  =  5  were  suggested  by  Wilson  (1989). 

Based  on  experimental  evidence,  Asano  (1992)  suggested  that 
the  erosion  depth  dt  (defined  as  the  difference  between  the  initial 
bed  level  and  the  still  bed  level  after  sheet  flow  started)  is  governed 
by  the  dynamic  Coulomb  criterion,  which  states  that  a  shear  stress 
and  its  normal  stress  are  proportional.  The  normal  stress  at  the 
bottom  consists  of  static  pressure  of  the  particle  lattice,  dispersive 
pressure  due  to  particle  collisions,  and  pore-fluid  turbulence  stress. 
Applying  this  criterion  at  the  boundary  between  the  mobile  and 
immobile  layers  yields  the  following  relationship  (the  bed  shear 
stress  is  defined  at  the  still  bed  level  contrary  to  Dick  &  Sleath), 
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fO 

'w.  6„  =  /  flg(s-  l)Ctan0rdz  (24) 

J—db 

where  q>r  is  the  critical  dynamic  angle  of  internal  friction.  Assuming 
that  C  and  (pr  are  constant  over  the  layer,  the  erosion  depth  can  be 
related  to  the  Shields  parameter  6: 


d$o  Ctan  (pr 

Approximate  estimates  of  C  and  (pr,  C-0.4  and  <pr  =  30°  in  Eq.  (25), 
yields: 

=  5  e  (26) 

U50 

Since  6b  <  <5S,  and  as  Wilson  (1989)  deduced  as  from  ao,  this  means 
that  the  coefficient  as>  1,  which  yields  a  value  that  is  at  least  twice 
as  large  as  the  one  proposed  by  Wilson  (1989).  It  should  be  noted 
that,  depending  on  the  way  the  sheet-flow  layer  is  estimated  (visual 
observations,  concentration  profile  including  the  “tail”,  or  linearly 
extrapolated  to  zero),  strong  differences  may  be  observed.  Indeed, 
Sumer  et  al.  (1996)  estimated  the  thickness  of  the  sheet-flow  layer 
from  visual  observations  and  from  the  concentration  profiles.  Using 
the  second  method,  values  of  ds  were  twice  as  large  as  values 
obtained  using  the  first  one.  The  value  as  =  \  is  used  for  this  study. 

3.1.5.  Bed-load  sediment  transport 

Another  method  to  estimate  the  roughness  height  is  proposed 
here.  If  the  net  sediment  transport  rate  is  known  (over  a  half-cycle 
or  a  full  cycle)  and  assuming  a  quasi-steady  behavior  for  the  sedi¬ 
ment  transport  (no  phase  lag  between  the  fluid  velocity  and  the 
sediment  concentration),  the  effective  Shields  parameter  can  be 
derived  employing  a  sediment  transport  formula  that  describes  the 
relationship  between  the  Shields  parameter  and  the  sediment 
transport  rate  over  flat  beds.  It  should  be  noted  that  this  method  to 
estimate  the  roughness  height  would  be  biased  if  strong  wave 
acceleration  or  boundary  layer  streaming  is  present.  These  effects 
may  indeed  influence  the  roughness  height  (Nielsen,  2002;  Myr- 
haug  and  Holmedal,  2005).  Nielsen  (1992)  proposed  to  use  the 
Meyer-Peter  and  Muller  (1948)  formula,  but  in  the  present  study, 
the  Soulsby  (1997)  formula  was  preferred  since  it  was  based  on  the 
quasi-steady  integration  of  the  Nielsen  formula  over  a  half  wave 
cycle  (c/.  Camenen  and  Larson,  2005), 

4s, 1/2  =  aw\Jg(s  -  l)d503(£-  dcr)3/2  (27) 

where  aw  =  5.1,  6  is  the  mean  value  over  a  wave  half-period  of  the 
instantaneous  Shields  parameter  (B  =  0.5 6  for  sinusoidal  waves), 
and  dcr  is  the  critical  Shields  parameter  for  incipient  movement  of 
the  sediment.  Even  tough  both  the  Meyer-Peter  &  Muller  and 
Soulsby  formulas  originally  included  the  grain  roughness  only, 
Camenen  and  Larson  (2005)  showed  that  these  formulas  per¬ 
formed  well  for  sheet-flow  transport  in  steady  flows  using  the  total 
or  equivalent  bed  shear  stress.  Once  the  Shields  parameter  has  been 
determined  from  the  sediment  transport  data,  the  friction  factor 
and  the  roughness  height  can  easily  be  estimated  using  Eq.  (10). 

For  the  full  cycle  data  (asymmetric  waves),  the  net  measured 
sediment  transport  is  obtained  from  qs ,  net  =  Qs,  1/2,  crest  -Qs,  1/2, 
trough*  Where  q5t  1/2,  crest  and  qs,  1/2,  trough  are  calculated  from  Eq.  (27) 
using  #1/2, crest  and  0i/2, trough*  respectively.  The  time-averaged 
Shields  parameters  are  based  on  the  assumption  that  the  friction 
factor  is  constant  over  the  wave  period  (discussed  in  Camenen  and 
Larson,  2005).  Thus,  the  parameter  values  are  only  functions  of  the 
velocity  time  series.  Assuming  a  wave  velocity  profile  following 
Stokes  2nd-order  wave  theory  (i.e.  uw(t)  =  Uw[cos  cot  +  rw  cos  2 cot], 


where  rw  is  the  wave  asymmetry  and  w  =  2tc/Tw)  with  a  half-period 
equal  to  Tw/2  for  both  crest  and  trough,  the  following  relationship  is 
obtained  (see  appendix  in  Camenen  and  Larson,  2007b): 


01/2, crest  _  01/2, trough  /^qn 

1  1  13 r  T  r2  ~  1  _  13 r  ir2~U 
1  Hh  37r  w  +  rw  1  37r'w  +  rw 

After  simple  algebraic  manipulations,  qS}  net*  to  the  first  order  of 
Ocr/B  and  rw,  is  approximated  by: 

qs, net ~—rw(iw^g(s  -  l)d5O303//2  (  29) 

Eq.  (29)  is  valid  only  if  the  bed  shear  stress  largely  exceeds  its 
critical  value  for  the  inception  of  transport  and  rw  <  1 . 


3.2.  Data  sets 

Table  2  summarizes  the  compiled  data  sets,  where  the  type  of 
flow  motion  (experimental  set-up),  the  methods  used  to  estimate 
the  roughness,  the  number  of  data  points,  the  sediment  properties 
(material  used,  density,  median  grain  size)  as  well  as  the  range  of 
values  for  the  main  hydrodynamic  parameters  (wave  orbital 
velocity  and  period,  skin  Shields  parameter)  are  listed. 

It  can  be  observed  that  most  of  the  data  are  from  Oscillating 
Water  Tunnels  (OWT).  This  kind  of  experiment  has  one  advantage 
for  this  study:  large  orbital  velocities  and  Shields  parameters  can  be 
reached  (sheet- flow  regime).  Previously,  experimental  studies 
were  often  carried  out  using  an  Oscillating  Tray  (OT;  oscillating  bed 
in  a  tank  of  still  water,  cf.  Sleath,  1978).  More  recently,  some 
experiments  by  Dohmen-Janssen  and  Hanes  (2002)  were  carried 
out  in  a  Large  Wave  Flume  (LWF).  Also,  light  material  was  often 
used  for  experiments  on  sheet-flow  regime  and  velocity  profiles, 
whereas  sand  was  generally  used  for  experiments  on  sediment 
transport. 

All  data  where  ripples  were  observed  were  discarded.  For  many 
cases  where  the  Shields  parameter  is  relatively  low,  ripples  are 
expected  to  occur.  However,  for  all  the  case  used  in  this  study,  the 
ripples  did  not  have  time  to  develop  (half-cycle  experiments).  In 
case  of  the  data  set  by  Watanabe  and  Isobe  (1990),  the  observed 
plane  bed  may  be  because  the  measurements  were  carried  out 
before  the  ripples  develop. 

In  Eq.  (27),  it  is  assumed  that  the  contribution  from  suspended 
load  amounts  to  only  a  small  (i.e.,  negligible)  fraction  of  the  total 
load.  Based  on  experiments,  Dohmen-Janssen  (1999)  estimated 
that  suspended  load  typically  consists  of  approximately  20%  of  the 
total  load  for  high-stress  data  in  OWT  experiments,  which  is  in 
agreement  with  this  assumption.  It  should  also  be  noted  that  for 
the  method  using  the  bed-load  transport  formula  over  a  full  cycle, 
phase-lag  effects  can  occur  and  influence  the  result.  The  Dohmen- 
Janssen  (1999)  criterion  has  been  used  in  order  to  eliminate  the 
data  for  which  significant  phase-lag  occurs,  where  a  steady  equa¬ 
tion  (Eq.  (27))  for  the  bed  load  is  not  valid, 


dsoj 

Pp‘  =  w:>05 


(30) 


where  ppi  is  a  phase-lag  parameter  and  Ws  is  the  sediment  fall 
speed. 


3.3.  Validation  of  analysis  methods 

As  mentioned  above,  the  “velocity  profile”  method  may  be 
considered  as  the  “reference  method”.  Thus,  using  data  sets  where 
the  “velocity  profile”  method  and  other  method  can  be  employed,  it 
is  possible  to  make  comparisons  for  validation  of  the  methods. 
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Table  2 

Data  summary  for  oscillatory  flow  under  plane-bed  conditions  (OWT:  Oscillating  Water  Tunnel;  OT:  Oscillating  Tray,  LWF:  Large  Wave  Flume). 


3.3 A.  Momentum  equation  in  the  fluid/ sediment  mixture 

Dick  and  Sleath  (1991)  estimated  the  bed  shear  stress  using  both 
“velocity  profile”  and  “momentum  equation”  methods.  A  comparison 
between  the  two  methods  (see  Fig.  la)  produced  very  good  agree¬ 
ment  after  the  shear  stress  estimated  from  the  “momentum  equation” 
method  was  multiplied  by  2.  The  underestimation  due  to  the  choice  of 
initial  bed  level  for  calculating  the  bed  shear  stress  or  the  influence  of 
sand  concentration  on  the  flow  only  included  in  the  “momentum 
equation”  method  may  explain  this  systematic  difference. 

3.3.2.  Energy  dissipation 

For  the  “energy  dissipation  factor”  method  proposed  by  Nielsen 
(1992),  no  other  measurements  were  found  to  validate  this  method. 
Wijetunge  and  Sleath  (1998)  measured  both  friction  and  energy 
dissipation  coefficients  at  the  levels  z  =  0  and  z  =  -ft.  The  friction 
coefficient  was  estimated  from  the  momentum  equation  in  the 
fluid/sediment  mixture.  It  appears  that  fe  and  fw  are  somehow 
related;  but  the  assumption  fe  ~fw  is  not  validated. 

3.3.3.  Sheet-flow  layer  thickness 

Several  authors  found  a  relationship  between  the  erosion  depth 
db  and  the  Shields  parameter  based  on  theoretical  or  experimental 
investigations  (see  Asano’s  relationship,  Eqs.  (25)  and  (26)), 

db  =  ab0^bd5O  (31) 

with  (ab,  /3b )  =  (3,  3/4)  following  Sawamoto  and  Yamashita  (1986), 
(ab,  ft>)  =  (8.5,  1)  following  Asano  (1992),  and  (cq,,  ft,)  =  ( 3,  1) 
following  Zala  Flores  and  Sleath  (1998).  Dohmen-Janssen  (1999) 
proposed  different  values  depending  on  the  grain  size,  i.e.,  (at, 
ft, )  =  (7.8, 1 ),  if  d50  <  0.2  mm,  and  {ab,  fib)  =  (3.5, 1 ),  if  d$ o  >  0.2  mm. 
Flowever,  most  of  these  authors  computed  the  Shields  parameter 
using  the  skin  friction. 


Fig.  2a  depicts  ft/dso  versus  the  total  Shields  parameter,  where 
the  roughness  height  is  estimated  as  ks  =  ft.  Although  the  uncer¬ 
tainties  remain  quite  large,  most  of  the  data  points  are  found  to  be 
surrounded  by  the  equations  proposed  by  Asano  and  Zala  Flores  & 
Sleath  (The  Asano  equation  tends  to  overestimate  the  results  as  he 
used  the  skin  Shields  parameter  to  fit  it,  which  leads  to  a  horizontal 
shift).  As  observed  by  O’Donoghue  and  Wright  (2004),  the  corre¬ 
lation  between  the  erosion  depth  and  the  Shields  parameter  seems 
to  occur  after  a  critical  value  for  the  Shields  parameter  corre¬ 
sponding  to  the  inception  of  sheet  flow  ( 6cr ,  s/~  0-8).  Using  this 
methodology,  it  seems  that  the  erosion  depth  (and  the  sheet-flow 
layer  thickness)  is  proportional  to  the  Shields  parameter  to  the 
power  fib ,  with  0.75  <  fib  <  1,  except  for  the  Wijetunge  &  Sleath  data 
set  where  fib  =  1.5.  This  means  that  the  previous  hypothesis  ( ks  =  db) 
is  reasonable. 

In  Fig.  2b  is  5sldb  plotted  versus  the  total  Shields  parameter 
where  the  roughness  is  estimated  as  ks  =  ds.  It  can  be  observed 
that  the  sheet-flow  layer  thickness  tends  to  be  equal  to  the 
erosion  thickness  when  the  Shields  parameter  reaches  a  critical 
value  varying  from  1  to  10.  The  four  points  (from  Sawamoto  & 
Yamashita  data  set)  where  the  ft/ft-value  is  large  correspond  to 
cases  where  suspension  occurred,  which  markedly  affected  the 
results.  The  data  from  Dick  &  Sleath,  where  ds  <  db,  can  also  be 
explained  since  the  sand  may  have  been  arranged  loosely  before 
the  experiment  and  packing  occurred  when  the  sheet-flow 
regime  appeared.  For  lower  shear  stress,  a  maximum  value  is 
reached  of  ft/ ft,  =  3-4.  It  may  correspond  to  an  intermediate 
regime  at  the  beginning  of  the  sheet-flow  regime,  where  insta¬ 
bility  effects  on  the  movement  of  each  particle  are  larger  (vertical 
movement  not  negligible),  and  where  the  bed  level  slightly  rises. 
Assuming  suspended  sediment  is  negligible,  by  conservation  of 
sediment,  the  relation  ft /ft  =  Cmax/C  is  obtained.  Assuming 
C  =  0.4  and  Umax  —  0.6,  d  s/db  =  1.5  would  be  expected.  Based  on 
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Fig.  1.  Validation  of  the  “momentum  equation”  method  (a)  and  “sheet-flow  thickness”  method  (b)  against  the  “velocity  profile”  method  (solid  line  corresponds  to  a  perfect 
agreement,  dashed  lines  to  an  agreement  within  a  factor  2;  data  from  Dick  and  Sleath,  1991  in  (a);  DS:  Dick  and  Sleath,  1991,  WS:  Wijetunge  and  Sleath,  1998,  DJ:  Dohmen-Janssen 


et  al.,  2001  in  (b)). 


this  assumption  and  Fig.  2,  the  sheet-flow  layer  seems  to  be  loose 
(C~  0.1— 0.2)  for  relatively  low  shear  stress  and  is  getting  very 
dense  (C~Cmax)  for  large  shear  stress. 

As  a  first  approximation,  it  may  be  assumed  that  ks  =  5S.  It  may 
be  noted  that  Sumer  et  al.  (1996)  observed  a  similar  relationship  as 
Eq.  (31 )  for  steady  flow  with  (cq,,  ft,)  =  (5, 1 ).  For  the  same  data  set, 
the  equation  proposed  by  Wilson  (1966)  that  ks  =  56d 50,  shows 
good  agreement.  This  means  ks  ~  5S 

Dick  and  Sleath  (1991),  Wijetunge  and  Sleath  (1998)  and  Doh¬ 
men-Janssen  et  al.  (2001)  estimate  the  bed  shear  stress  using  the 


“velocity  profile”  method  and  measured  the  sheet-flow  layer 
thickness  (or  erosion  depth).  A  comparison  between  the  results 
obtained  employing  the  two  methods  (see  Fig.  lb)  presents  a  larger 
discrepancy  than  for  the  “momentum  equation”  method.  An 
overestimation  of  the  results  is  often  observed  mainly  for  the  data 
from  Dick  &  Sleath.  However,  this  could  be  explained  similarly  as 
for  Fig.  2a;  the  Dick  &  Sleath  data  set  corresponds  to  the  highest 
values  of  ft,/d5 0  observed.  It  should  be  remembered  that  some 
discrepancy  may  result  from  the  way  the  sheet-flow  layer  thickness 
is  estimated. 
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Fig.  2.  (a)  Erosion  thickness  5 b  and  (b)  the  ratio  <5 s/<5b  versus  the  total  Shields  parameter  calculated  using  the  relationship  ks  =  8b. 


3.3.4.  Bed-load  sediment  transport 

The  validation  of  the  “bed-load  transport”  method  is  not  as  easy  as 
for  the  previous  methods  since  no  data  set  was  found  where  both 
sediment  transport  under  waves  only  and  the  velocity  profiles  were 
measured.  Nnadi  and  Wilson  (1992)  performed  experiments 
for  steady  currents  in  pressurized  closed  conduits  where  the  rough¬ 
ness  height  ks  could  be  computed  using  both  the  classical  Nikuradse 
equation  (assuming  a  logarithmic  velocity  profile)  and  employing 
measured  bed-load  transport  rate.  The  Meyer-Peter  &  Muller  formula 
was  utilized  to  estimate  the  roughness  height  from  these  bed-load 
data.  As  observed  by  Bayram  et  al.  (2003)  (see  Fig.  1),  these 
two  different  methods  are  highly  correlated,  especially  for  the 
experimental  data  with  sand.  However,  the  “bed-load  transport” 
method  seemed  to  yield  results  that  depend  on  the  sediment  density. 
Through  this  validation,  it  is  however  believed  that  the  two  methods 
used  here  to  estimate  ks  under  oscillatory  flow  are  compatible. 

Finally,  the  data  sets  from  Sawamoto  and  Yamashita  (1986), 
Asano  (1995),  Dohmen-Janssen  and  Hanes  (2002),  Ahmed  and  Sato 
(2003),  and  O’Donoghue  and  Wright  (2004)  present  measurements 
of  both  the  sheet-flow  layer  thickness  and  the  net  sediment  trans¬ 
port.  Assuming  that  the  “sheet-flow  layer”  method  gives  accurate 
results  with  ds  =  ks,  it  is  possible  to  validate  the  “bed-load  transport” 
method.  Unfortunately,  for  the  data  sets  of  Ahmed  &  Sato  and 
O’Donoghue  &  Wright,  since  strong  phase-lag  effects  occurred,  only 
a  limited  portion  of  the  data  could  be  used  since  the  Soulsby  formula 
(1997,  cf.  Eq.  27)  for  the  estimation  of  the  bed  load  was  not  valid 
anymore.  The  selection  of  the  data  points  was  made  using  the 
relationship  by  Dohmen-Janssen  (1999)  that  yields  2izdsl 
( WsTw )  >  0.35  (limit  above  which  quasi-steady  models  are  not  valid). 

Fig.  3  illustrates  variations  of  ks  and  6  using  both  “bed-load 
transport”  and  “sheet-flow  layer  thickness”  methods.  Dohmen- 
Janssen  and  Hanes  (2002)  showed  that  the  thickness  of  the  sheet- 
flow  layer  varies  with  time  and  measured  ds  over  a  wave  period.  The 
<5s-values  observed  at  the  crest  and  at  the  trough  of  the  waves  were 
plotted.  The  induced  Shields  parameters  correspond  to  the 
maximum  value  at  the  crest  and  at  the  trough,  defined  as  follows, 

n  _  V^/w, crest/trough uw, crest/trough2 

crest/trough  -  (S-l)gd50  K  ’ 

where  fWt  Crest  and  fw<  trough  are  the  friction  coefficients  calculated 
using  Eq.  (10)  with  ks  =  5s,  crest  and  ks  =  5s,  trough  respectively.  The 
results  observed  in  Fig.  3  show  fairly  good  agreement  between 
the  two  methods,  though  with  marked  dispersion.  However,  for  the 
Sawamoto  &  Yamashita  data  set,  it  seems  that  a  lower  value  on  at, 


would  imply  better  results.  This  can  easily  be  explained  as  Ss/3t,  ~  3 
for  most  of  the  cases  in  this  data  set  (see  Fig.  2b).  Moreover,  the 
large  underestimation  observed  for  the  remaining  Ahmed  &  Sato 
and  O’Donoghue  &  Wright  data  may  be  caused  by  the  phase-lag 
effects  occurring  in  these  data.  Indeed,  the  criterion  proposed  by 
Dohmen-Jansen  (cf.  Eq.  (30))  is  nearly  fulfilled  (0.3  <ppi<  0.5)  for 
most  of  these  cases. 

4.  Roughness  height  for  plane  bed  under  oscillatory  flow 

4.1.  Comparison  between  existing  roughness  height  relationships 
and  data 

In  order  to  make  the  results  clearer  and  to  understand  the  limits 
of  each  method  used,  the  data  have  been  plotted  in  three  separate 
figures  depending  on  the  method  employed.  In  Fig.  4a  are  the 
results  from  the  “velocity  profile”,  “momentum  equation”  and 
“energy  dissipation”  methods  plotted;  in  Fig.  4b  from  the 
“boundary  layer”  method;  and  in  Fig.  4c  from  the  “bed-load 
transport”  method.  The  three  figures  present  the  effective  rough¬ 
ness  ratio  versus  the  total  Shields  parameter  based  on  the 
“measured”  ks.  Also,  the  semi-empirical  relationships  from  the 
literature  to  predict  roughness  height  are  shown  in  all  the  graphs. 
The  Grant  &  Madsen  equation  (1982)  could  not  be  plotted  on  the 
same  graph,  since  it  depends  on  the  skin  Shields  parameter. 

In  Fig.  4a,  it  is  clear  that  the  data  points  and  the  fitted  curve  from 
the  “energy  dissipation  method”  proposed  by  Nielsen  (1992, 
pp.152-155)  show  a  totally  different  behavior  than  the  other  data 
sets  and  formulas.  However,  the  larger  roughness  heights  obtained 
by  Nielsen  might  be  explained  by  the  fact  that  the  plane  regime  was 
unstable  (the  measurements  were  taken  before  ripples  had  time  to 
form).  As  shown  in  Section  3.1.2,  the  assumption  fe  ~/w  seems  not 
to  be  satisfactory  (at  least  for  this  case).  These  data  points  will  thus 
be  rejected  in  the  following  study. 

The  equation  proposed  Xu  and  Wright  (1993)  tends  to  over¬ 
estimate  the  roughness  height  for  most  of  the  cases.  The  Wilson 
(1989)  and  Ribberink  (1998)  formulas  yield  much  more  accurate 
results.  The  slope  (3  (see  Eq.  (1))  observed  from  the  data  points 
appears  to  be  higher  than  1  (/?  ~  1.5  using  a  least-square  fit). 
Furthermore,  the  data  from  Sleath  (1987)  confirm  that  a  minimum 
value  of  ks/dso  exists  when  the  Shields  parameter  is  lower  than  1 
(no  sheet  flow),  corresponding  to  a  fixed  bed.  However,  the  classical 
value  observed  by  Kamphuis  (1974)  and  Yalin  (1977),  and 
confirmed  by  Camenen  et  al.  (2006),  i.e.,  fcs/dso  =  2,  overestimates 
the  results  for  this  data  set. 
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Fig.  3.  Validation  of  the  “bed-load  transport”  method  against  the  “sheet-flow  thickness”  method  (solid  line  corresponds  to  a  perfect  agreement,  dashed  lines  to  an  agreement 
within  a  factor  2). 


The  roughness  heights  obtained  using  the  “sheet-flow  layer” 
method  (see  Fig.  4b)  display  a  slightly  different  behavior.  Indeed, 
the  observed  slope  is  (3  ~  1,  but  this  could  easily  be  explained  by  the 
basic  assumption  of  this  method,  i.e.,  ks  <5S  (Eq.  (23)),  since  <5S  (or  8b) 
is  proportional  to  the  Shields  parameter  to  the  power  1  (c/.  Eqs.  (25) 
and  (31)  and  Fig.  2).  It  should  also  be  noted  that  the  data  from 
Dohmen-Janssen  and  Elanes  (2002)  obtained  in  a  large  wave  flume 
show  larger  roughness  height  values  for  a  fixed  Shields  parameter 
compared  to  the  data  from  OWT.  Elowever,  Dohmen-Janssen  & 
Elanes  showed  that  the  experimental  set-up  does  not  really  influ¬ 
ence  the  thickness  of  the  sheet-flow  layer.  The  main  difference 
appeared  to  be  the  way  to  estimate  the  sheet-flow  layer  thickness: 
they  calculated  it  from  the  sediment  concentration  profiles,  while 
previous  authors  used  visual  estimations.  Moreover,  Ahmed  and 
Sato  (2003)  used  a  PIV  system  to  predict  the  sheet-flow  layer 
thickness  that  apparently  also  displayed  larger  observed  values  on 
the  sheet-flow  layer  thickness  for  a  fixed  Shields  parameter. 

Finally,  the  results  obtained  through  the  “bed-load  transport” 
method  (see  Fig.  4c)  tend  to  confirm  those  from  “velocity  profile” 
and  “momentum  equation”  methods:  the  slope  (3  determined  by 
least-square  fitting  toward  the  data  points  appears  to  be  higher 
that  1  (/?  ~  1.5).  Also,  the  observed  roughness  heights  are  gener¬ 
ally  larger  than  for  the  previous  methods  when  6  >  1.  When  6  <  1, 
the  obtained  results  are  very  scattered  with  ks/dso  varying  from 
0.1  to  10.  For  some  data  from  Watanabe  &  Isobe,  the  Shields 
parameter  being  very  close  to  its  critical  value,  the  results 
obtained  using  Eq.  (29)  may  be  also  distorted  (underestimated). 
Elowever,  Camenen  et  al.  (2006)  observed  a  similar  dispersion  for 
the  estimation  of  the  roughness  height  under  steady  currents  and 
for  relatively  low  shear  stress. 


oscillatory  sheet  flow:  “Mobile  beds  at  high  shear  stress  are  neither 
smooth  boundaries  nor  rough  ones,  but  obey  their  own  frictional 
law,  analogous  to  the  other  cases  but  with  length  scale  proportional 
to  the  sheet-flow  thickness  <5S.  It  follows  that  sheet-flow  behavior 
can  be  made  formally  equivalent  to  the  rough-boundary  case  by 
setting  the  equivalent  Nikuradse  sand-grain  roughness  ks  equal  to 
a  multiple  of  <5S”  for  both  steady  and  oscillatory  flows.  Ribberink 
(1998)  followed  this  idea  and  suggested  that  a  constant  represen¬ 
tative  value  for  ks  may  be  used  during  the  wave  cycle.  Ele  proposed 
a  similar  formula  for  steady  and  oscillatory  flows  (apart  from  the 
grain-related  or  skin  roughness  that  he  assumed  different  for 
a  steady  flow,  ks  =  3dgo,  and  for  an  oscillatory  flow,  ks  =  dso),  that  is 
function  of  the  Shields  parameter  and  the  sediment  grain  size.  Zala 
Flores  and  Sleath  (1998)  also  observed  that  the  flow  may  be 
considered  as  quasi-steady  for  the  sheet-flow  regime  at  relatively 
low  shear  stresses  as  pressure  gradient  and  inertia  forces  are 
negligible.  They  found  that  the  sheet-flow  thickness  may  thus  be 
proportional  to  the  Shields  parameter  times  the  medium  grain  size. 
When  the  instantaneous  shear  stress  reaches  its  maximum,  if  there 
is  no  phase-lag  between  the  flow  and  the  shear  stress,  the  accel¬ 
eration  is  zero  and  the  shear  stress  should  act  similarly  to  a  steady 
flow.  This  means  that  the  effective  roughness  ratio  may  be  a  func¬ 
tion  of  the  maximum  wave-induced  Shields  parameter  in  the  same 
way  as  for  a  steady  current. 

Thus,  following  the  results  from  Camenen  et  al.  (2006),  it  is 
hypothesized  that  the  relationship  could  be  written  /<s/d50  «  01,7  as 
soon  as  the  Shields  parameter  reaches  a  critical  value.  The  value 
(3  =  1.7  is  quite  close  to  the  slope  observed  for  the  data  in  the  graphs 
of  Fig.  4.  A  new  equation  for  the  roughness  height  for  plane  bed 
under  oscillatory  flow  is  proposed, 


4.2.  A  new  relationship 


-r-  =  a  +  b 

“50 


(33) 


4.2.1.  Analogy  to  steady  current 

Based  on  a  study  on  the  vertical  velocity  profile  above  a  sheet- 
flow  layer,  Wilson  (1989)  observed  an  analogy  between  steady  and 


where  a  and  b  are  coefficients  ((a,  b)  =  (0.6, 2.4)  for  steady  flows,  cf. 
Camenen  et  al.,  2006)  and  0cr>  ur  corresponds  to  the  limit  between 
the  lower-regime  and  the  upper-regime.  The  upper-regime  may  be 
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In  Fig.  5  (where  all  the  data  were  used  apart  from  the  “energy 
dissipation”  method  since  it  was  not  validated  against  any  other 
method  and  the  results  are  quite  different  from  the  other  methods), 
the  effective  roughness  ratio  is  plotted  versus  the  total  Shields 
parameter  with  the  different  sediment  densities  used  for  the 
experiments  emphasised.  It  may  be  noted  that  the  differences 
observed  in  Section  4.1  may  not  only  be  due  to  the  different 
methods  used  to  compute  the  roughness  height  but  also  due  to  the 
different  kind  of  sediments  used  in  the  various  experiments. 
Indeed,  it  appears  that  in  all  the  experiments  where  the  “velocity 
profile”  method  was  applied,  the  sediment  density  was  particularly 
low.  This  may  imply  that  the  effective  roughness  ratio  is  not  only 
a  function  of  the  Shields  parameter  but  also  a  function  of  the  sand 
characteristics  (such  as  the  sediment  density)  or  the  flow 
characteristics. 


4.2.2.  Critical  Shields  parameter  for  the  upper-regime 

The  main  objective  of  the  present  study  is  to  develop  a  new 
predictive  formula  for  the  roughness  height  that  can  perform 
satisfactory  under  a  wide  range  of  conditions.  A  sensitivity  analysis 
was  made  to  estimate  the  influence  of  the  main  parameters,  which 
are  the  median  grain  size  dso,  the  settling  velocity  Ws,  the  sediment 
density  s,  the  wave  period  Tw,  and  the  orbital  velocity  at  the  bottom 
Uw,  in  Eq.  (33).  In  order  to  keep  the  analogy  with  steady  flows,  an 
equivalent  Froude  number  for  waves  Fj-w  and  the  dimensionless 
settling  velocity  Ws*  proposed  by  Camenen  et  al.  (2006)  were 
introduced, 


Frw  — 


Uw 

Cg5 


(34) 


where  5  =  ^/vTw  is  the  thickness  of  the  viscous  (Stokes)  layer,  v  the 
kinematic  viscosity  of  the  fluid,  and: 


Ws*  = 


(5-1) 

(gv) 


2lV3 


ws 


(35) 
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Fig.  4.  Effective  roughness  ratio  versus  total  Shields  parameter,  together  with  semi- 
empirical  predictive  relationships,  including  data  obtained  using  (a)  the  “velocity 
profile”  (V.P.),  “momentum  equation”  (M.E.)  and  “energy  dissipation”  (E.D.)  methods; 
(b)  the  “sheet- flow  layer”  method;  and  (c)  the  “bed-load  transport”  method. 


defined  as  the  regime  where  ks  is  no  more  a  function  of  the  grain 
size  only.  On  the  other  hand,  the  lower-regime  may  be  defined  as 
the  regime  ks  is  only  a  function  of  the  grain  size  (skin  friction).  If 
6  <  Scr,  ur,  ks  <  3dso,  which  corresponds  approximately  to  the  skin 
friction. 


A  study  of  the  root-mean-square  error  using  Eq.  (33)  with  varying 
powers  on  Frw  and  Ws *  was  carried  out.  Finally,  by  introducing  the 
critical  Shields  parameter  for  upper  plane-bed  regime  6cr,  ur,  such  as 
/<s/d50  =  3  when  6  =  0cr>  ur  (in  order  to  obtain  an  average  of  ks/dso 


Fig.  5.  The  effective  roughness  ratio  versus  the  total  Shields  parameter  (using  all  data 
apart  from  the  “energy  dissipation”  method)  including  the  influence  of  sediment 
density  (the  dotted,  dashed,  solid,  and  dashed-dotted  lines  correspond  to  Eq.  (33)  with 
a  =  0.6  and  b  =  2.4,  and  0cr,  ur  =  0.25,  0.5, 1.0  and  2.0,  respectively). 
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when  8  <  0cr,  ur  approximately  equal  to  2),  the  best  fit  was  found  for 
a  —  0.6,  and  b  =  2.4  with: 


=  0.115- 


171.20 
1  rw 


W°*40(s  -  1 


\0.30 


(36) 


Apart  from  the  stronger  effect  of  the  relative  density,  Eq.  (36)  is 
quite  similar  to  the  equation  obtained  for  the  critical  Shields 
parameter  for  the  upper-regime  under  steady  flows  (i.e.,  8cr< 
ur=l .1 2F?40 Ws^070).  The  present  results  appear  to  be  not  as 
sensitive  as  the  steady  flow  results  to  the  settling  velocity  but  more 
sensitive  to  the  sediment  density. 


4.2.3.  Comparison  with  experimental  data 

Fig.  6  presents  a  comparison  between  effective  roughness  ratios 
ks/dso  obtained  from  the  data  and  corresponding  values  computed 
with  Eq.  (33).  It  appears  that  Eq.  (33)  is  not  markedly  affected  by 
the  different  methods  used  to  estimate  the  roughness  height. 
Nevertheless,  the  effective  roughness  ratio  was  found  to  be 
proportional  to  the  Shields  parameter  to  the  power  1.7,  whereas 
the  sheet-flow  layer  thickness  appeared  to  be  proportional  to  the 
Shields  parameter  to  the  power  0.75-1.5  (see  Section  3.3.3).  Thus, 
the  hypothesis  for  the  “sheet-flow  layer  thickness”  method  that  is 
ks  =  5s  seems  to  be  accurate  enough  even  if  it  corresponds  to  the 
maximum  value  for  the  proportionality  with  the  Shields  param¬ 
eter.  Indeed,  if  Eq.  (33)  is  used  to  predict  5S  (assuming  ks  =  8S)  more 
them  95%  of  the  data  are  correctly  predicted  within  a  factor  2. 
Some  data  points  derived  using  the  “sediment  transport”  method 
are  also  overestimated  by  Eq.  (33),  but  these  data  correspond  to 
Shields  parameter  values  close  to  the  critical  value  for  incipient 
motion  where  the  method  is  expected  to  induce  larger 
uncertainties. 

Table  3  summarizes  the  predictions  of  equivalent  roughness 
height  for  the  different  formulas  within  a  factor  of  2  (P2)  and  5  (P5) 
of  the  measured  values  (“factor  x”  means  between  x  times  and  1/x 
times  the  measured  roughness  height),  as  well  as  the  mean  value 
( mf )  and  standard  deviation  {sf)  of  the  ratio /=  ks<  pred IK  expe-  The 
data  set  used  for  these  calculations  corresponds  to  all  previously 
presented  data  apart  from  the  data  which  clearly  produce  inaccu¬ 
rate  results  (that  is  data  based  on  the  “energy  dissipation”  method 
and  the  data  from  the  “sediment  transport”  method  where  strong 
phase-lag  were  observed).  It  can  be  noted  that  there  is  a  poor 


Fig.  6.  Comparison  between  effective  roughness  ratio  obtained  from  the  data  and 
corresponding  values  computed  with  Eq.  (33)  (the  bold  line  corresponds  to  a  perfect 
agreement  and  the  dashed  lines  correspond  to  prediction  within  a  factor  of  2). 


Table  3 

Prediction  of  equivalent  roughness  within  a  factor  of  2  and  5  of  the  measured  values 
together  with  a  logarithmic  error  index  for  the  studied  formulas  and  all  data. 


agreement  between  all  previous  formulas  and  the  data  within 
a  factor  of  2  (less  than  50%  of  the  data),  and  they  sometimes  yield  an 
error  of  up  to  one  order  of  magnitude.  Thus,  a  logarithmic  error 
index  is  proposed  to  establish  a  quantitative  measure  of  the 
agreement  between  the  formulas  and  the  data.  The  term  “En0g” 
(logarithmic  error  index)  is  defined  as, 


ttiog  — 


_  100y^ 

~  n  4^ 


logio 


(/<s/d5o)pred 
{ks  /d$o)  expe 


(37) 


in  which  n  denotes  the  number  of  observations.  An  error  of  a  factor 
10  would  thus  lead  to  Er\og  =  100. 

Large  overestimations  are  observed  for  the  Grant  &  Madsen, 
Nielsen  and  Xu  &  Wright  formulas.  However,  it  should  be  noted  that 
the  data  set  from  Nielsen  was  not  used  for  this  comparison,  which 
explains  the  poor  results  for  this  formula.  The  Wilson  and  Ribberink 
formulas  tend  to  slightly  underestimate  the  results  for  large  Shields 
parameters  (0  >  1).  However,  for  6  <  1,  the  Wilson  formula  under¬ 
estimates  the  observation  since  it  yields  no  minimum  value.  The 
Ribberink  formula  also  produces  underestimations  since  its 
minimum  value  ( ks  =  dso)  is  generally  smaller  than  the  observed 
values.  Eq.  (33)  produces  predictions  that  are  much  more  accurate 
than  previous  relationships  for  most  of  the  data  sets.  From  Table  3, 
it  is  clear  that  Eq.  (33)  gives  the  best  results  for  all  measures  of 
accuracy  employed.  Eq.  (33)  was  however  calibrated  using  the  same 
data  set  it  was  tested  against. 

Fig.  7  displays  the  significant  improvement  of  the  roughness 
height  predictions  using  Eq.  (33)  (bold  lines)  compared  to  formulas 
depending  on  the  Shields  parameter  only  (the  Wilson  formula  is 
added  as  a  reference,  dashed  lines).  In  order  to  construct  these 
graphs,  the  data  were  divided  into  8  groups  for  0  <  Ocr,  ur  <  0.20, 
0.20  <  dcr,  ur  <  0-40,  0.40  <  8cr,  ur  <  0-60,  0.60  <  8cr,  ur  <  0.80, 
0.80  <  8cr,  ur  <  1.00,  1.00  <  8cr,  ur  <  1-50,  1.50  <  8cr,  ur  <  2.00,  and 
2.00  <  8cr,  ur  <  2.50.  Equation  (33)  was  plotted  using  the  mean  value 
of  8cr,  ur  for  each  group. 


4.2.4.  Inception  of  sheet-flow  and  upper-regime 

The  inception  of  sheet  flow  corresponds  to  a  situation  where  the 
wave  ripples  are  disappearing,  simultaneously  as  the  energy  is 
increasing  (increasing  wave  orbital  velocity).  Thus,  a  relationship 
should  exist  between  the  inception  of  sheet  flow  and  the  inception 
of  the  upper-regime  defined  in  the  previous  section.  Camenen  and 
Larson  (2006)  proposed  a  maximum  value  of  the  wave  orbital 
velocity  Uw<  Crsf  for  the  inception  of  the  sheet- flow  regime  or  wash¬ 
out  of  the  wave  ripples, 

Uw.crsf  =  8.35^J  (s  —  V)g\J  d505w(l  +  rw)  (38) 

where  rw  =  (uw,  max  +  Ww,  min)/(2Uw)  is  the  wave  asymmetry.  A 
critical  Shields  parameter  8cr,  sf  for  the  inception  of  sheet-flow  may 
be  obtained  from  Eq.  (38)  adding  a  coefficient  which  is  a  function  of 
the  dimensionless  grain  size  d*  =  ^/(s  -  I)g/v23d50  (c/.  Hanson 
and  Camenen,  2007): 
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Fig.  7.  Dimensionless  roughness  height  versus  Shields  parameter  for  different  values  on  the  critical  Shields  parameter  for  the  upper-regime  (the  bold  lines  correspond  to  Eq.  (33) 
and  the  dashed  lines  to  the  Wilson  formula,  Eq.  (3)). 


6crj  =  10d;3/4,fel  +rw)2  (39) 

y  °w 

The  Van  Rijn  formula  (1993)  which  is  only  a  function  of  d*  ( 6cr , 
Sf=  26 dcr)  yields  very  similar  results  though  it  is  not  a  function  of 
the  wave  period. 

In  Fig.  8  are  the  different  equations  for  the  critical  Shields 
parameters  for  the  inception  of  sheet  flow  and  upper-regime 
plotted  versus  the  dimensionless  grain  size.  As  8cr<  ur  is  a  function 
of  the  wave  orbital  velocity,  it  has  been  plotted  for  four  different 
values  of  Uw  from  0.25  to  2  m/s  (with  Tw  =  6s  and  rw  =  0).  It 
appears  that  0cr,  ur  yields  a  similar  behavior  as  6cr,  sf  with  respect 
to  the  grain  size  (decreasing  function  of  d*).  Eqs.  (36)  and  (39) 


Fig.  8.  Comparison  between  the  critical  Shields  parameter  for  the  upper-regime  and 
critical  Shields  parameters  for  the  inception  of  sheet  flow  (VR:  Van  Rijn,  1993;  C&L: 
Camenen  and  Larson,  2005,  2006,  2007a, b). 


also  have  a  similar  sensitivity  to  the  wave  period  (or  <$w). 
However,  Eq.  (36)  for  6cri  ur  differs  from  the  two  other  relation¬ 
ships  for  6cr,  sf  as  it  is  also  sensitive  to  the  wave  orbital  velocity 
and  then  yields  much  smaller  value  when  Uw  <  0.5  m/s.  By  defi¬ 
nition,  the  lower  plane-bed  regime,  where  the  equivalent 
roughness  height  may  be  approximate  as  the  grain-related 
roughness  height  (i.e.  ks  ~  2dso),  is  differentiated  from  the  upper 
plane-bed  regime  for  which  ks  >  2dso.  Thus  the  sheet-flow 
regime  is  a  specific  case  of  the  upper  plane-bed  regime  (for  the 
sheet-flow  regime,  ks  was  always  observed  to  be  much  larger 
than  2d5o)  and  0cr<  ur  <  6Cr,  sf 

5.  Conclusion 

In  the  present  study,  data  from  various  sources  were  compiled 
and  analyzed  to  determine  the  roughness  height  for  plane-bed 
regime  under  oscillatory  flow.  Based  on  this  analysis,  the  following 
conclusions  were  drawn: 

•  Several  methods  were  used  and  proposed  to  estimate  the 
roughness  height.  The  “velocity  profile”,  “momentum  equa¬ 
tion”  and  “sediment  transport”  methods  appear  to  give  a  fairly 
consistent  estimation  of  ks.  On  the  other  hand,  the  “energy 
dissipation”  method  induces  larger  values  on  the  roughness 
height.  However,  this  method  was  not  validated  against  any 
other  method. 

•  The  sheet-flow  layer  thickness  was  found  to  be  a  function  of 
the  total  Shields  parameter.  The  ratio  between  the  sheet-flow 
layer  thickness  and  erosion  depth  varies  from  3  for  medium 
shear  stresses  to  1  for  large  shear  stresses.  The  assumption  that 
ks  =  ds  appeared  to  be  correct. 

•  Flow  resistance  for  plane  beds  under  oscillatory  flow  can  be 
expressed  in  terms  of  the  ratio  of  Nikuradse’s  equivalent  sand 
roughness  to  the  grain  diameter,  fe/dso,  given  as  function  of  the 
Shields  parameter  based  on  the  total  shear  stress,  and  a  critical 
Shields  parameter,  according  to  Eq.  (33).  0cr,  ur  defines  the  limit 
between  the  lower-regime  where  ks/d50  may  be  consider  as 
constant  and  the  upper-regime  where  fcs/dso  is  a  function  of  the 
Shields  parameter. 
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•  Similar  results  for  the  effective  roughness  height  were 
observed  for  the  oscillatory  flow  as  for  the  steady  current  (c/. 
Camenen  et  al.,  2006),  i.e.,  /<s/d5o  is  proportional  to  the  total 
Shields  parameter  to  the  power  1.7  as  soon  as  6cr,  ur  is  reached. 

•  The  critical  Shields  parameter  for  the  upper-regime  0cr,  ur  was 
found  to  be  a  function  of  the  settling  velocity,  the  sediment 
density,  and  a  wave-equivalent  Froude  number.  This  value  is 
strongly  related  to  the  critical  Shields  parameter  for  inception 
of  sheet  flow,  but  may  yield  also  much  lower  values.  As 
a  result,  it  appears  that  the  roughness  height  may  increase  as 
soon  as  sediment  transport  occurs.  Also,  the  sheet-flow 
regime  was  found  to  be  a  specific  case  of  the  upper  plane-bed 
regime. 

•  The  proposed  relationship  (Eqs.  (33)  and  (36))  yields  the  best 
results  among  the  studied  formulas.  Also  these  equations  may 
also  be  used  to  estimate  the  sheet-flow  layer  thickness  as  soon 
as  the  sheet-flow  regime  is  reached. 

•  A  relationship  between  the  roughness  height  and  the  total 
Shields  parameter  makes  iterative  computations  necessary  in 
practical  applications.  Such  a  method  could  imply  large 
discrepancies  in  the  results  compared  to  the  direct  fitting  using 
measured  data,  and  thus  needs  to  be  improved. 

•  Some  other  effects  may  influence  the  roughness  height  and 
friction  coefficient.  Myrhaug  and  Holmedal  (2005)  argued  that 
the  boundary  layer  streaming  may  affect  the  bed  shear  stress 
for  laminar  (with  relatively  high  Reynolds  numbers)  and 
smooth  turbulent  flows  (with  relatively  low  Reynolds 
numbers).  Acceleration  does  affect  the  instantaneous  shear 
stress  (and  roughness  height)  by  inducing  a  phase  lead  and  an 
increase/decrease  of  the  shear  stress  for  accelerating/deceler¬ 
ating  flows.  Thus,  the  instantaneous  roughness  height  (and 
friction  coefficient)  may  fluctuate  around  its  mean  value. 
Although  these  fluctuations  may  induce  some  onshore- 
directed  sediment  transport  (Nielsen,  2002;  Camenen  and 
Larson,  2007a),  the  purpose  of  this  study  was  to  better  identify 
the  mean  value  on  the  roughness. 
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